Numerical study of the topological Anderson insulator in HgTe/CdTe quantum wells 
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We study the disorder effect on the transport properties in the HgTe/CdTe semiconductor quan- 
tum wells. We confirm that at a moderate disorder strength, the initially un-quantized two terminal 
conductance becomes quantized, and the system makes a transition to the novel topological An- 
derson insulator (TAI). Conductances calculated for the stripe and cylinder samples reveal the 
topological feature of TAI and supports the idea that the helical edge states may cause the anoma- 
lous quantized plateaus. The influence of disorder is studied by calculating the distributions of local 
currents. Base on the above-mentioned picture, the phenomena induced by disorder in the quantum 
spin Hall region and TAI region are directly explained. Our study of the local current configurations 
shed further light on the mechanism of the anomalous plateau. 
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I. INTRODUCTION 

Quantum spin Hall effect (QSHE), which is proposed 
as a new class of topological state of matter in two dimen- 
sions, has generated a great deal of interest 1 . Contrary 
to the integer or the fractional quantum Hall state, which 
is induced by a magnetic field that breaks the time re- 
versal symmetry, QSHE is caused by a strong spin-orbit 
interaction that maintains the time reversal symmetry. 
Kane and Mele proposed a Z2 classification for this kind 
of the new topological stated. For materials with spatial 
inversion symmetry, the index could be easily computed 
as the product of parities of the wave function at sev- 
eral high symmetry points in the Brillouin zone. The Z2 
classification can be generalized to three-dimensional sys- 
tems as well^. Recently, topological insulators suitable 
for room temperature applications are also predicted 4 . 
The QSH state has the helical edge states, namely, hav- 
ing two counter-propagating edge states for the two op- 
posite spin polarizations. The helical edge states are sta- 
ble against time-reversal conserving perturbations, since 
backscattering processes need to connect the upper and 
lower edges of the sample. The probability of backscat- 
tering is exponentially suppressed as the sample width is 
increased. Recent experiment^ provides evidences of the 
QSHE in HgTe/CdTe quantum well (QW) structures, as 
predicted theoretically 6 . The de-coherence effect in QSH 
samples is also investigated^ Some interesting questions 
emerge such as the response of a QSH state to the disor- 
der, and the interplay of the helical edge states and the 
bulk states. 

According to the celebrated scaling argument, in two 
dimensions the bulk electron wave functions are local- 
ized in the presence of any weak disorder- . But there 
are two known exceptions, one is the two dimensional 
systems with a strong spin-orbital coupling, and the 
other is the quantum Hall transition between different 
plateaus^. With the discover of QSHE state in two di- 
mensions, study of the localization is in demand. Sheng 



et al investigated the disorder effect in honeycomb lat- 
tice QSHE sample, they found the QSHE phase is robust 
against weak disorder 10 . Ono and Obuse et al studied 
the critical properties of the transition from a metal to 
a QSH state, and they found the results to be somewhat 
controversia l 11 ' 12 . 

Recently, Lii£ et al studied the transport properties of 
the HgTe/CdTe QWs in the presence of disorder, they 
found as the increasing of the disorder strength, the ini- 
tially un-quantized conductance became quantized, i.e. 
the sample enters into the QSH state because of the dis- 
order, so they named the state as "topological Anderson 
insulator" (TAI). The newly anomalous quantized con- 
ductance plateau is caused by the edge transport, which 
is indirectly revealed by the unchanged plateau value 
along with width variation for two-terminal calculation 
and quantized transmission coefficient for four-terminal 
calculation. However, the detailed mechanics of the edge 
transport is less clear. 

In the present paper, we study the effect of disorder on 
the electronic state of HgTe/CdTe QWs. We carry out 
Keldysh's nonequilibrium-Green-function (NEGF) calcu- 
lations based on a four-band tight binding model. First, 
we perform the conductance calculations for two differ- 
ent geometries. In the case of a stripe geometry(see fig 
.1(a)), the presence of the helical edge states are evident 
from the band spectrum, the TAI phase described in ref- 
erence 12 appears. While for the cylindrical geometry, 
i.e. periodical boundary condition along y direction(see 
fig. 1(b)), there is no edge state, the bulk state is local- 
ized by disorder and there is no quantized conductance. 
These results strongly support the thesis that anomalous 
conductance plateau is due to the edge transport and 
give a better understanding of the topological feature 
of the TAI phase. Second, we obtain the distributions 
of the local currents for the two-terminal strip samples 
with different chemical potentials and disorder strengths. 
The evolvement of the local current vector configurations 
gives rise to a direct demonstration of the impurity influ- 
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ence in the HgTe/CdTe QWs. By analyzing these local 
current configurations, the transport phenomena in both 
normal QSH region and TAI region are clearly explained. 
Moreover, the detailed results, such as the coexistence of 
the bulk and edge states at the dip point (see Fig. 6) and 
the bulk states assisted backscattering obtained from the 
local current vector configurations shed further light on 
the mechanism of the disorder induced edge states for the 
TAI. 

The rest of this paper is organized as follows. In Sec- 
tion II, we introduce the effective tight-binding model. 
The formulas and calculation method are also described. 
In Section III, the numerical results and their discussions 
are presented. Finally, a conclusion is given in Section 
IV. 



II. THE MODEL AND METHOD 

As a starting point, we introduce the effective Hamil- 
tonian for the HgTe/CdTe QWs with Anderson impurity 
in the tight-binding representation. We consider a square 
lattice with four special orbit states |s, t),bx + ipy, T) , 
I s ; l)> I — (Px — iPy)-, i> on eacn site. Here t,J, denotes 
the electron spin. Through symmetry consideration, the 
effective Hamiltonian can be written as£ & 13 ' 14 
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Here i = {ix, iy) is the site index, and Sx and Sy are unit 
vectors along the x and y directions, (fi = (ctj, a, b\, di) T 
represents the four annihilation operators of electron on 
the site i with the state indexes |s, t)JPx + ipy, T) , I s , I): 
— (Px — iPy), l>) respectively. The on-site matrix ele- 
ments satisfy E isa — E s + W isa , E ipc = E p + W ipc , E is b = 
E s + Wisb-,Eip d = E p + W ipd . Wisa, W ipb , W^c and 
Wipd are on-site disorder energies uniformly distributed 
in the range [— with the disorder strength W. 

E s ,E p , V S s, Vpp, and V sp are the five independent param- 
eters that characterize the clean HgTe/CdTe samples. It 
is clear that near the T point the lattice Hamiltonian (Eq. 
(1)) in k-representation can be reduced to the continu- 
ous Hamiltonian in Ref.6 when we take V sp = —iA/2a, 
Vss = (B + D)/a 2 , Vpp = (D - B)/a 2 , E s — C + M — 



4(B + D)/a 2 , and E p = C — M — A(D — B)/a 2 . Here a is 
the lattice constant and all the parameters A, B,C, D, M 
can be controlled experimentally 5 . Moreover, the Eq.(l) 
can be directly obtained by discretizing spatial coordi- 
nates of the continuous Hamiltonian using the substitu- 
tion k x > — igj, ky > — i^. 15 

In this paper, we apply the model to two geometric de- 
vices. The device a (see fig. 1(a)) is of stripe geometry, 
while the device b (see fig. 1(b)) is of cylindrical geome- 
try which can be obtained by rolling the device a into a 
tube. Without disorder(W=0), the energy spectrum of 
such two geometric devices can be calculated by diago- 
nalizing Eq.(l) using the periodic boundary condition in 
x directio n 16 ' 17 . Next we investigate how the disorder af- 
fects the transport properties of such systems. For both 
devices a and b, the size of the central region is L x W. 
For convenience, we assume that the Anderson impurities 
only exist in the (red) filled region and the temperature 
is zero. 

In our simulations, a small external bias V = Vl — Vr is 
applied between the two terminals. With the help of the 
NEGF method, the local current between neighboring 
sites i and j is calculated from the formul a 18 ' 19 

Ji^j = ^-ME H ^Gf^ ia (E F )](V L - V R ) (2) 

a, 13 

where Vl,Vr are the voltages at the Lead-L,R. a, (3 de- 
note the state indices. G n (E F ) = G r T L G a is electron 
correlation function with line width function T a (Ep) = 
i[Tr a {E F ) - (T, r a {E F ))+] , the Green function G r (E F ) = 
[E F - H cen - Yj r L {E F ) - Yj^Ep)]- 1 and the Hamilto- 
nian in the central region denoted as H cen . The retarded 
self-energy E£,S^j due to the coupling to the lead-L, R 
can be calculated numerically 20 . Note that for Eq. (1), 
the spin up(ai,Ci) subsystem and spin down {b\,di) sub- 
system are decoupled. The local current between neigh- 
boring sites i and j with spin index a can also be 
calculated from Eq.(2) by summing over only the state 
index with the corresponding unitary subsystem. The 
current Jj, flowing through the device is calculated by 
summing all the local current Ji—,i + sx for an arbitrary 
column. After obtaining the current Jl , the linear con- 
ductance G LR is given by G L r = Jl/(Vl - Vr). In ad- 
dition, the linear conductance can be directly obtained 
by G LR = Tr[T L G r T R G a }. The agreement between the 
two methods gives strong confirmation of our analytical 
derivations and numerical calculations. 

In the following numerical calculations, we choose the 
realistic material parameters A, B, C, D, M that arrived 
from the HgTe/CdTe QWs.- The sample width (or cir- 
cumference) is fixed to L y = 80a with the lattice constant 
a = 5nm . Since the model is only valid in small fc, we 
confine the Fermi energy within a small region near the 
r point. In the presence of disorder, the conductance G, 
conductance fluctuation 8G, the local current etc., 
are all averaged over up to 500 random configurations. 
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FIG. 1: (color online) (a) and (b)are the schematic diagrams 
for two devices. The Anderson impurities only exist in the 
(red) filled regions, (c) and (e) show the one-dimensional 
energy bands for device a. The parameters are A = 3.645eV ■ 
A,B = -68.6eV • A 2 ,C = 0.0meV,D = -51.2eV • A 2 , and 
gap parameter M = —lOmeV (c), M — 2meV (e). (d) and 
(f) demonstrate the one-dimensional energy bands for device 
b with the same parameters as (c) and (e), respectively. 
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FIG. 2: (color online) The conductance G (a)(b) and con- 
ductance fluctuation SG (c)(d) vs disorder strength W for 
different Fermi energy Ef and gap parameter M. (a)(c) for 
device a and (b)(d) for device b. The central region length 
L x = 200a. Other parameters (A, B, C, and B) are the same 
as fig. l(c)-l(f). 



III. NUMERICAL RESULTS 

Let us first compare the two devices in geometry and 
topology. In device a with open boundary (fig. 1(a)), 
there are two edges in the y direction. Thus, the edge 
states can exist in such device. Since the fig. 1(b) is a 
cylinder without edge, the edge states are prohibited in 
the device b. In figures 1(c) and 1(e), the band struc- 
tures of clean HgTe/CdTe QWs of are plotted. There 
exist a bulk energy gap approximately of 2\M\ in both 
figures. Moreover, there are two degenerate bands (edge 
states) cross inside the gap for gap parameter M < (fig. 
1(c)). In contrast, the crossing bands vanished when M 



was tuned up to above zero (fig. 1(d)). These results 
are in agreement with those of previous studies^. In 
figures 1(d) and 1(f), we plot the corresponding band 
structures in the cylindrical geometry. For both gap pa- 
rameter \M\ > and \M\ < 0, though the degeneracy is 
higher due to the enhanced symmetry, the bulk energy 
gap is nearly unchanged and there is only tiny shift in 
the energy bands. However, there is a big difference be- 
tween the two samples, the edge states which cross inside 
the gap vanish in the later one. The phenomenon origi- 
nates from the topology of the device, when the sample 
is sufficient large, the bulk states are hardly affected by 
the topology change at the edge while the edge states 
are totally destroyed. In conclusion, rolling the sample 
from a strip to a cylinder destroys the edge channels but 
maintain the bulk state properties. 

Next, we investigate how the transport properties are 
affected by disorder. Fig .2(a) and 2(c) plot the conduc- 
tance G and SG versus disorder strength for device a. 
When the system is in the inverted regime (M < 0) with 
the Fermi energy inside the bulk gap (Ep = 7meV), for 
a range of disorder strength W, the two terminal conduc- 
tance is quantized without much fluctuations. Such ob- 
servation agrees with the previous result that the QSHE 
is robust against weak disorder^iS. However, when the 
chemical potential is tuned up into the bulk band region 
near the T point, no matter M > or M < 0, the new 
intriguing phenomena emerge. The conductance G de- 
creases while the fluctuation SG increases when the dis- 
order is first applied. When the disorder strength contin- 
ues to increase, instead of localization, the conductance 
begins to increase to a quantized value (2e 2 /h) and main- 
tains at this value for a certain range before eventually 
decreases. Meanwhile, the conductance fluctuation SG 
decreases to zero and vanishes for the corresponding W. 
The anomalous conductance plateau indicates that the 
sample becomes a topological insulator. More impor- 
tantly, different from the traditional topological insula- 
tor (QHE etc), the quantized value is induced by the im- 
purities. This disorder induced anomalous conductance 
plateau was discovered in a very recent workJ^ Natu- 
rally, there exists a question: what is the mechanism that 
causes such an anomalous plateau? 

To answer this question, we first examine G and SG 
as a function of disorder strength for device b using 
the same parameters as device a [see Figs. 2(b) and 
2(d)]. Unlike the device a, the transport properties fol- 
low the traditional Anderson metal-insulator transition 
feature^ and the anomalous conductance plateau is ab- 
sent. For instance, the conductance monotonously de- 
creases to zero with increase of the disorder strength 
and the critic disorder strength W c increases with raising 
of the Fermi energy. Significantly, the metal-insulator 
transition point for device b is roughly at the start- 
ing point of the anomalous plateau for device a. Take 
AI = — lOmeV, Ep — l8meV for example [see black solid 
line in fig. 2(a) and 2(b)], the anomalous plateau sets 
up at the disorder strength W c ~ 93meV for device a 
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FIG. 3: (color online) Configurations of the 
flow vector for device a with Fermi energy Ef z 
lOmeV, central region size L x = 200a, L y - 
disorder strength W=0 (a) W=50meV(b) W= 
W=220meV (d). The inset in fig. 3(a) is the 
local current flow vector. The vector direction 
local current flow direction and the vector len£ 
tional to the logarithm of local current value. 
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and for that threshold value the device b becomes an 
insulator(G ~ 0,SG ~ 0). As shown in fig .1, the bulk 
states in the device a and the device b are the same. 
Only the edge states are completely destroyed by rolling 
the device a into the device b. This gives a direct evi- 
dence that the anomalous conductance plateau originates 
from the edge states. 

To get a better insight into the microscopic origin of 
the conductance variations, we examine the disorder ef- 
fect through the local-current-vector-flow configurations. 
Due to the time reversal symmetry, we only consider the 
spin up subsystem, the influence of spin down subsys- 
tem can be directly obtained by time reversal symmetry. 
Here, the local current flow vector on site i is defined as 
Jl = J} 



+i+ix 



J i^i+ay 



In fig. 3, the typical distributions of local currents for 
device a in traditional QSHE region are plotted. For 
a clean sample (see fig. 3(a)), the local currents locate 
mainly on the upper edge and their values decay expo- 
nentially towards the bulk. Surprisingly, the small dis- 
order initially makes the edge channel narrower (see fig. 



80 

60 

40 

20 


80 

60 

40 

20 


80 

60 

40 

20 


80 

60 

40 

20 





:4.(a) 



:4.(b): 



r<f\ ) ■ . 



4.(c): 



■ i , / v \ \ >^ > -* \v 1 — -»* t - 



-\r> ' ^- 



4.(d) 



(0 



40 



80 120 

X/a 



160 



200 



FIG. 4: (color online) Configurations of the local current flow 
vector for device a with the same sample sizes as for fig.3, pos- 
itive gap parameter M = 2meV , Fermi energy Ef = 18meV , 
and disorder strength W=0 (a) W=100meV (b) W=150meV 
(c) W=250meV (d). 



3(b)). Though the mechanism is unclear, we note such 
phenomenon was also observed recently by Chu et al£^ 
In their paper, this narrowing effect was indirectly ob- 
served by the decreasing of the oscillation period of the 
A-B ring, while in this paper such effect is directly shown 
by the spatial distributions of the local currents. When 
the Anderson disorder strength is getting larger, the lo- 
cal currents spread to the bulk and broaden the edge 
channels again(see as fig. 3(c)). However, only when 
the disorder strength exceeds the critical value W c , the 
spread local-current flow can reach the lower edge chan- 
nels with different chirality, the effective backscattering 
(as shown in the local current flow vector located near the 
lower edge in the region < X < 40a in fig .3(d)) can 
take place, leading to the reduction of the conductance 
between the two terminals. These pictures explain why 
the traditional quantized plateau is robust under weak 
disorder and how it is destroyed in the strong disorder 
limit. 

Next, the Fermi energy is tuned to Ep = 18meV, sit- 
ting slightly above the bulk gap. A positive gap parame- 
ter M = 2meV, for which there is no helical edge states 
inside the bulk gap for the clean sample, is chosen in 
the following simulations. The conductance G vs. dis- 
order strength W shown in fig. 2(a) can be classified to 
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FIG. 5: (color online) Configurations of the local current 
flow vector for device a with the same sample sizes as for 
fig .3, negative gap parameter M — — lOmeV, Fermi energy 
Ef = 18meV, and disorder strength W=0 (a) W=65meV (b) 
W=100meV (c) W=250meV (d). 
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FIG. 6: (color online) (a) Distribution of local currents in the 
disorder region for device a with length L x — 400a and L y = 
80a at the dip point, (b) The position-related current J p vs 
longitude axis x. J p is defined as the summation of jj_ i+Sx for 
four layers in the corresponding region. The other parameters 
are the same as (a), (c) The two-terminal conductance G vs. 
disorder strength W with different sample lengths L x . 



four regions (i) without disorder, (ii) before the anoma- 
lous plateau, (iii) on the anomalous plateau, and (iv) 
after the anomalous plateau. The typical configurations 
of local-currcnt-flow vector in such four regions are plot- 
ted in fig. 4. For the clean sample (see fig. 4(a)), the 
local current not only flow forward along the upper edge, 
but also uniformly flow forward in the bulk. While W 
is increased from zero into the region (ii), the local cur- 
rents in the bulk of the disorder region become smaller 
and more irregular in directions (see fig. 4(b)), which 
directly shows the decline of the bulk transport and the 
localization of the bulk states. Fig .4(c) demonstrates 
the distribution of the local currents when the device 
is in the region (iii). The most interesting phenomena 
exist in the disorder region X € [40a, 160a]. The local 
currents in the bulk decline to zero while the residual 
currents flow in the upper edge and flow without much 
scattering. In addition, throughout the region (iii), the 
bulk transport vanishes due to the localization of all the 
bulk states. In contrast, the edge transport shows the 
same behavior as in the traditional QSHE region: the 
local currents deeply spread into the bulk without any 
effective backscattering with W increasing until it con- 
nected the opposite edge channels(see as fig. 4(d)). Ob- 
viously, such edge transport will lead to the quantized 
plateau. In the negative gap parameter situation, for 
example M = — lOmeV, Ep = 18me^,(see fig. 5) the 
evolvement of local current configurations resembles that 
of M > situation. So far, the fig. 4(c) and fig. 5(c) 
are the most strong evidence directly showing that the 
TAI is caused by the edge transport. In addition, these 
plots provide a vivid microscopic picture demonstrating 
the influence of disorder on TAI. 

Up to now, we have explained the majority phenomena 
emerged in the TAI. Nevertheless, the transport proper- 
ties in the region (ii) of G — W figure demands a detailed 
study for the following two reasons, (i) The conductance 
is not monotonously decreasing with increasing disorder 
strength but shows a dip feature prior to the anoma- 
lous plateau. Obviously, at the dip point, the system is 
neither in TAI phase nor in normal Anderson insulator 
phase due to its nonzero and non-quantized conductance 
value, (ii) For both types of gap parameter M > and 
M < 0, the dip does exist prior to the anomalous plateau 
and the conductance G behaves similarly after the dip. 
Thus, revealing the cause for the dip feature may help 
us to understand the mechanism of the formation of the 
edge states. 

In fig. 6 (a), the distribution of the local currents with 
the disorder strength W being fixed at the dip is plotted. 
The local currents flow are larger close to the edges than 
in the bulk. The predominant edge transport is clearly 
seen. In fig. 6(b) we plot the position-related current 
J p versus longitude axis x. Where J p is defined as the 

summation of jl_^- l+Sx f° r f° ur l a y ers m the correspond- 
ing region. For example, given sample width L y = 80a, 

J p for the upper edge is defined as Xay=7? 3i->i+8x m With 
the help of J pi one can quantitatively analyze the local 
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currents. On the whole, the behavior of the position- 
related current J p is similar to the local current config- 
urations, but it is smoother with the position x. More 
significantly, one can observe from fig. 6(b), the local 
currents flow from left to right with rapidly decreasing 
magnitude for upper edge, vice versa for the lower edge. 
This phenomenon can be attributed to the bulk state as- 
sisted backscattering between two edges. The fact that 
the position-related current J p in the bulk is small but 
nonzero indicates that the bulk states are not fully lo- 
calized at the dip. The scattered carriers in the upper 
edge can hop through such delocalizcd bulk states to the 
lower edge which leads to the backscattering processes. 
In other words, disorder not only destroys the bulk trans- 
port but also quickly destroys the edge transport. Be- 
cause of this, the conductance G is lower than plateau 
value at the dip point. For the region after dip (see fig. 
2(a)), increasing of W slowly destroys the bulk trans- 
port, but it also suppresses the bulk-assisted backscat- 
tering mechanism at the same time. It leads to an over- 
all increase of the conductance G. When the bulk states 
become fully localized, the anomalous plateau shows up. 
In fig. 6(c), the conductance G versus W with different 
sample lengths L x are plotted. The dip feature is clearer 
for a large L x because of the increasing probability of the 
backscattering between the two edges. However, for all 
L x , the anomalous plateaus appear with the same disor- 
der strength due to the fully localized bulk states. 

IV. CONCLUSIONS 

In summary, the disorder effect in HgTe/CdTe quan- 
tum wells is studied. We confirm the existence of the 



topological Anderson insulator (TAI) phase. Conduc- 
tances calculated for the stripe and cylinder samples re- 
veal the topological feature of TAI and support the idea 
that helical edge states cause the anomalous quantized 
plateau. With the help of local-current-vector config- 
urations for different chemical potentials and disorder 
strengths, the basic physical phenomena emerged in the 
normal QSHE region and in the TAI region are clearly 
understood. In particular, the analysis of the local cur- 
rent configurations provides us with the importance of 
the bulk-states-assisted backscattering in TAI that in 
turn help us to understand the mechanism of the for- 
mation of the disorder induced edge states. 
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